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Originally designed as a new nuclear reactor monitoring device, the Nucifer detector has suc¬ 
cessfully detected its first neutrinos. We provide the second shortest baseline measurement of the 
reactor neutrino flux. The detection of electron antineutrinos emitted in the decay chains of the 
fission products, combined with reactor core simulations, provides a new tool to assess both the 
thermal power and the fissile content of the whole nuclear core and could be used by the Inter¬ 
national Agency for Atomic Energy (IAEA) to enhance the Safeguards of civil nuclear reactors. 

Deployed at only 7.2 m away from the compact Osiris research reactor core (70 MW) operating 
at the Saclay research centre of the French Alternative Energies and Atomic Energy Commission 
(CEA), the experiment also exhibits a well-suited configuration to search for a new short baseline 
oscillation. We report the first results of the Nucifer experiment, describing the performances of the 
~ 0.85 m® detector remotely operating at a shallow depth equivalent to ~ 12 m of water and under 
intense background radiation conditions. Based on 145 (106) days of data with reactor ON (OFF), 
leading to the detection of an estimated 40 760 7'e, the mean number of detected antineutrinos is 
281 ± 7(stat) ± 18(syst) Pe/day, in agreement with the prediction 277 ± 23lze/day. Due to the large 
background no conclusive results on the existence of light sterile neutrinos could be derived, how¬ 
ever. As a first societal application we quantify how antineutrinos could be used for the Plutonium 
Management and Disposition Agreement. 


I. INTRODUCTION 

In a context of increasing needs for carbon emission- 
free energy, civilian nuclear power generation has the po¬ 
tential to play an important role in global energy pro¬ 
duction, and the list of countries aiming to acquire tech¬ 
nological know-how in the field of civilian nuclear en¬ 
ergy could increase. As a consequence, the International 
Atomic Energy Agency (IAEA) has been evaluating the 
potential of new technologies to guarantee that nations 
use nuclear energy only for peaceful purposes. 

Neutrino detectors have the unique ability to non- 
intrusively monitor a nuclear reactor’s operational status, 
thermal power and fissile content in real-time, from out¬ 
side the reactor containment. More specifically the sce¬ 
narios of interest are to confirm the absence of unrecorded 
production of fissile materials in declared reactors as well 
as to estimate the total burn-up[T] of a reactor core. Nu¬ 
cifer is a detector first built for long term reliable safe¬ 
guards measurements in the vicinity of operating nuclear 
reactor cores. The experiment aims at demonstrating 
the concept of “neutrinometry” at the pre-industrialized 
stage. Therefore a well-established technology and com¬ 
mercial components were chosen for the detection system. 


In addition, after the re-evaluation of antineutrino 
spectra mu, the reanalysis of all short baseline reac¬ 
tor experiments [5] lead to what is known as the Re¬ 
actor Antineutrino Anomaly (RAA) with = 

0.936 ± 0.024 [B]. And it turns out that the Nucifer ex¬ 
perimental set-up is adequate to probe the Reactor An¬ 
tineutrino Anomaly [5] , searching for possible oscillations 
into sterile neutrino species at very short baselines. In¬ 
deed, Nucifer is a compact detector (^ 1.2 m in diameter, 
0.7m in height), deployed at only 7.21 ± 0.11 m (centre to 
centre) from a compact nuclear core (57 x 57 x 60 cm®). 
Therefore the three conditions to search for short base¬ 
line oscillations are met. Nucifer is even the second short¬ 
est baseline reactor neutrino experiment, the only clos¬ 
est experiment being a 3.2L Gd-loaded liquid scintilla¬ 
tor detector deployed at only 6.5 m from the Savannah 
River Plant (USA) in 1965 [3 [5]. The third shortest 
baseline is the ILL experiment [9] with 377 L of liquid 
scintillator at 8.76 m from the core, which lead to a ra¬ 
tio of experimental to expected integral positron yield of 
0.832 ± 3.5 %(stat) ± 8.87 %(syst) as stated in the 1995 
reanalysis [TO] . 

Among the possible societal applications the Plu¬ 
tonium Management and Disposition Agreement 
(PMDA) [TT] could be monitored through neutrino 
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rate monitoring. Indeed, in this procedure, weapon- 
grade plutonium could be processed into Mixed OXide 
uranium-plutonium (MOX) fuel, irradiated in civil 
nuclear power reactors, and therefore transformed into 
material unusable in the fabrication of nuclear weapons. 

The paper first describes the experimental lay-out in 
part |nj from the Osiris reactor to the data acquisition 
system through the Nucifer detector. We then calculate 
the expected signal and its associated errors in part [ml 
accounting for the major part of the final uncertainty due 
to the relatively low detector efficiency. Then the cali¬ 
bration system is depicted in part EYi followed by the 
data analysis: we detail the data sample in section [V A| 
and the neutrino candidate selection in section V B then 
the accidental background in section [VCl the cosmic in¬ 


duced correlated background in section |V^ 


induced correlated backgrounds in section V 


the reactor 


^and finally 


the detection efficiency and the associated uncertainties 
in section m We thus present the Nucifer results at 
Osiris in part m summarized by the ratio of observed 
to expected neutrino detection rate. We finally discuss 
a potential application for the PMDA agreement, using 
Nucifer’s collected data, in part |VII[ 


II. EXPERIMENTAL LAYOUT 
A. Deployment site: the Osiris reactor 

Nucifer is installed on the concrete foundation slab of 
the Osiris reactor building, 11 meters beneath the wa¬ 
ter pool level, in a dedicated room next to the reactor 
core (see Fig. [^. Such a location allows to safely sup¬ 
port the weight of the detector together with its heavy 
passive shielding (> 62tonnes). It also offers a modest 
overburden, reducing the muon flux by a factor of 2.7 
with respect to sea level. The overburden is equivalent 
to ^ 12 meters of water. 

The reactor core, approximately at the Nucifer room 
ceiling level, is located 7.21 ± 0.11m away from the de¬ 
tector, from the centre of the core to the centre of the 
detector. In this configuration radiations from the re¬ 
actor core are attenuated by about 2 m of concrete and 
3.5 m of water. Nevertheless, the level of gamma radia¬ 
tions in the Nucifer room is still quite high during the 
reactor operations, enhancing the challenge of properly 
extracting the neutrino flux. 

Osiris is a light water experimental reactor of open- 
core pool type located at the Saclay research centre 
of the French Alternative Energies and Atomic Energy 
Commission (CEA). It operates at a nominal thermal 
power of 70 MW and has been designed for technologi¬ 
cal irradiation purposes and radioisotope production m- 
The fissile elements produce a high neutron flux in the 
core, at the level of a few 10^"^ neutrons cm“^s“^, both 
in the thermal and fast energy range. The core size is 
57 X 57 X 60 cm^, excluding its vessel, corresponding to 
56 cells loaded with 38 standard fuel elements (U 3 Si 2 Al 



FIG. 1. The Nucifer experimental layout. The detector centre 
is located 7.21 ± 0.11m away from the reactor core centre. 
East is on the right, south points to the reader. 


plates enriched at 19.75% in 6 control elements 

(made of Hafnium absorber in the upper part and of fuel 
in the lower part), 7 Beryllium elements used as neutron 
reflectors, and 5 cells equipped with water boxes dedi¬ 
cated to experiments. The cells are arranged in a square 
lattice with a lattice parameter of 8.74 cm. The absence 
of pressurization vessel allows for a direct access to the 
core at any time. 



FIG. 2. Top scheme of the Osiris core. “BG” squares (red) 
stand for core control elements, bottom core line squares 
(labelled 10, green) for Beryllium reflector elements. The 
square surrounding the core (purple) represents the steel ves¬ 
sel. Other areas (yellow, grey) stand for locations used for 
experimental irradiation devices, or industrial irradiation lo¬ 
cations m- North is on top, Nucifer is located on the left. 
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Osiris typically operates 180 days per year with 3 week 
reactor cycles, core being refuelled by about 1/7*^. The 
19.75% enrichment of the nuclear fuel in ^35^ 
short cycle duration suppress the evolution of the iso¬ 
topic composition of the nuclear fuel. To start the reac¬ 
tor, control rods 1, 2, 5 and 6 are first raised, and the 
core reaches criticality during the raising of the next rod 
which can be the 3 or the 4 (see fig . This rod is then 
raised to the top within the next few days, and the last 
rod is progressively raised during the remaining part of 
the cycle. This sequence leads to a displacement of the 
fission barycentre, in the east/west direction and verti¬ 
cally, different for the two rod configurations (3 then 4 
or 4 then 3). 


B. Nucifer detector 


The Nucifer detector (Fig. consists of a cylindrical 
liquid scintillator tank, surrounded by an active plastic 
scintillator veto to tag cosmic ray muons, and two lay¬ 
ers of shielding: 14 cm of boron-doped polyethylene to 
capture neutrons, and 10 cm of lead layer to attenuate 
external gamma rays. 

Three additional lead walls have been erected to fur¬ 
ther attenuate the reactor induced gamma rays. A first 
10 cm thick wall was originally installed between the de¬ 
tector and the reactor wall, on the east direction side 
(see Fig. [^. Then a second 4 cm thick wall was later 
built between the original lead wall and the reactor wall 
to further reduce the amount of gamma rays originating 
directly from the core (see section V CI resulting in a 
factor of 3 additional attenuation. 

Finally, a third 10 cm thick wall was erected on the 
southern side of the detector (see Fig. to suppress the 
gamma radiations coming from the primary water loop 
of the reactor cooling circuit located behind the 1 m thick 
concrete wall of the detector casemate. Indeed when flow¬ 
ing through the fuel elements, the water of this circuit is 
being highly activated by the fast neutron flux inside the 
core. is produced through (n,p) reaction on of 
the water m and subsequently decay with a 7.13 s half- 
life emitting a 6.1 MeV y-ray with an intensity of 67% 
(plus a 7.1 MeV y-ray with intensity 4.9 %). A dedicated 
circuit, located partly behind the southern wall of the Nu¬ 
cifer casemate, delays its arrival to the primary pumps 
by about 90s, giving time for ^®N to decay completely. 
The effect of this added shielding was to reduce the single 
event rate by more than an order of magnitude. 

The liquid scintillator is contained inside a stainless 
steel cylindrical vessel (1404.2 mm in height, 1250.4 mm 
in diameter). The internal surface of this vessel is 
coated with reflective white Teflon, chosen for its chem¬ 
ical compatibility with the liquid scintillator as well as 
for increasing visible light collection. The tank contains 
846.8 ± 7.0 L of liquid scintillator doped with a Gadolin¬ 
ium complex (a Gd-beta-diketonate) to enhance the cap¬ 
ture of thermal neutrons and sign the neutron capture. 


The chemistry of the Nucifer scintillator is based on the 
target liquid of the Double Ghooz experiment [15], with 
three main modifications. First the Gd-concentration 
was increased to 0.17% in mass to reduce the capture 
time of thermal neutrons (t„), expected to be in the range 
of 20 ps. Second the concentration of o-PXE (ortho- 
Phenylxylylethane) in the scintillator was increased to 
57 % in volume, the remaining 43 % being dodecane. 
In this way light yield and pulse shape discrimination 
power for rejection of correlated background events are 
improved. Finally, the concentration of the primary 
fluor PPO (2,5-Diphenyloxazole) was slightly increased 
to compensate for light losses due to the higher Gd con¬ 
centration m- However, the concentration of the sec¬ 
ondary fluor bis-MSB (l,4-Bis(2-MethylStyryl)Benzene) 
was kept at 20 mg/1. 

The light collection is performed by sixteen 8-inch 
photomultipliers (PMTs), type Hammamatsu R5912, lo¬ 
cated at the top of the detector vessel. A 25 cm thick 
acrylic disk separates the PMTs from the target liq¬ 
uid scintillator. Filled with mineral oil, this so-called 
buffer optically couples the PMTs to the liquid scintilla¬ 
tor. This design allows a good uniformity of the detector 
response to energy deposition in the whole target vol¬ 
ume and shields the scintillator from the intrinsic PMT 
radioactivity. The buffer is hxed to the tank lid. No 
phenomenon of light emission from PMTs was detected 
(so-called “flashing PMTs”). 

A light injection system guiding light from LED to 
Teflon diffusers in the tank through optical fibres allows 
the monitoring of the PMT gain as well as the liquid op¬ 
tical properties, thanks to several types of light diffusers. 
In addition small encapsulated radioactive sources can 
be deployed along the target central axis inside a vertical 
stainless steel tube externally coated with Teflon. 



(£) Stainless steel tank 
Teflon coated 

(D Nitrogen atmosphere 

® 16 PMTs 

® Acrylic buffer 

d) Target Liquid (0.8 m^) 

d) Calibration tube 

( 2 ) Cosmic ray veto 

® 14 cm B-doped 
polyethylene 

® 10 cm Lead 


FIG. 3. Cut view of the Nucifer detector. The overall volume 
is about 3 X 3 X 2.4 m®. 

The muon veto is made of 32 modular detectors, 
each one containing a 5 cm thick plastic scintillator bar 
(150 to 170 cm length, 25 cm width) associated to a single 
PMT decoupled from its surface. The plastic scintilla¬ 
tor thickness and light collection were optimized in order 
to discriminate cosmic muons from high energy gamma 
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rays in a compact and cost-effective way. Indeed, muon 
energy deposition is at least lOMeV, corresponding to 
perpendicularly crossing muons at the minimum ionizing 
point, while the highest gamma energy is also lOMeV, 
from neutron capture on metals surrounding the detector 
or in the concrete. Using events that saturate all PMTs 
in the Nucifer vessel and checking if the muon veto had 
triggered or not, we measured a muon detection efficiency 
of 97%. 


C. Electronics and data acquisition 

The analogical output of each of the 16 PMTs is split 
in 4 different channels. The first channel is used to build 
the trigger (the analogical sum of all PMTs), and the sec¬ 
ond channel is routed to a constant fraction discriminator 
to get the time information. The two last channels are 
delayed and sent to commercial CAEN Charge to Digital 
Converter (QDC) modules. One integration gate is ad¬ 
justed to integrate the full signal centred around a 150 ns 
time window. The other gate is delayed by about 40 ns 
to integrate only the late component of the signal. The 
comparison of the resulting two charges can be later used 
for Pulse Shape Discrimination (PSD). 

The acquisition is triggered either by the analogical 
sum of all PMTs overcoming a threshold equivalent to 
about 1 MeV, or by any of the muon veto module over¬ 
coming a threshold equivalent to about 10 MeV, or by 
computer driven LED and random signals. The dead 
time is computed online by counting internal clock pulses 
inhibited by all hardware vetos (such as QDC busy gates 
or saturation of the buffers of the QDC/TDC). 

The Data Acquisition system (DAQ) is based on the 
Lab VIEW software allowing a remote control of the ac¬ 
quisition and a constant monitoring of the safety pa¬ 
rameters, such as pressure, liquid level and temperature 
at various locations. The data stream is automatically 
chopped up in 500 MB runs and transferred to the CC- 
IN2P3 computing centre for storage and off-line analysis. 
Typical run duration is 43 minutes with the Osiris reactor 
operating at full power, and 53 minutes when the reactor 
is not operating. The difference comes from the reactor 
induced gamma background in the Nucifer tank. 

In case of a significant deviation of any set of predefined 
parameters, warning emails are automatically sent to on- 
call experts. If the deviation exceeds a higher threshold, 
an alarm is sent to the reactor control room. Since its 
installation at Osiris in Spring 2012, the Nucifer detector 
has been operating without any safety failure. 


liquid scintillator free protons: Ag -I- p —t e"*" -I- n, taking 
advantage of the time delayed coincidence between the 
positron and neutron signals. The positron detection is 
the prompt signal, its energy is related to the neutrino’s 


by: E. 


prompt 


E^ — 0.782 MeV. Then the neutron is 


captured with high efficiency either on Gd or H atoms 
in the liquid scintillator. Neutron captures occur with 
a mean time of about 20 ps and lead to the emission of 
a few y-rays with a total energy of 8.05 MeV in average 
for Gd natural isotopic composition and thermal neutron 
capture at ambient temperature. 

To first order, the number of neutrino events detected 
per day depends on the neutrino flux emitted by the reac¬ 
tor, the baseline, the number of target free protons, and 
the detector efficiency. The neutrino flux depends on the 
reactor thermal power and, since different fissioning iso¬ 
topes lead to different neutrino spectra, on the fission 
fractions (see Fig. 15). Finally, the event rate at a given 
time t is given by 


P (f) 

T^{t) = ’ - PpY^akit) a,-AE)Sk{E)kB dE 

Y,au{t)E^ V 

k 

( 1 ) 

The first term describes the number of fissions per unit 
time with Pth(t) the thermal power, ak{t) the fission 
fraction of isotope k {k =235U, 239p^^ 24ip^) 

and Ek the mean energy per fission of isotope k. Pp 
is the proton density in the liquid scintillator. The in¬ 
tegral on the antineutrino energy represents the mean 
fission cross-section with (Jcr^{E) the IBD cross section 
and Sk{E) the fission spectrum of isotope k in units of 
Vg fission" ^ MeV" ^ [IB] . 

The effective length is homogeneous to a distance 
and takes into account the finite extensions of both the 
core and the detector, a necessary refinement due to the 
very short baseline of the experiment. It is defined as 


^eff 



clMd^Tt 


with 'Lf,fc(t, ^) the fission density of isotope k, normal¬ 
ized to one hssion, Yg the reactor core volume, fd the 
detector volume, e{E,rd) the detection efficiency, rt a 
vector pointing to a given point in the core and ra a 
vector pointing to a given point in the detector. In the 
case of point-like core and detector one readily recovers 
the equivalence Pplgg = J^p/dirL'^ with Afp the number of 
protons in the target and L the mean baseline. 


III. EXPECTED SIGNAL 

A. Theoretical neutrino rate 

The detection of Osiris antineutrinos is achieved 
through the Inverse Beta Decay (IBD) reaction on the 


B. Parameters and associated uncertainties 

All relevant reactor operating parameters are being 
provided by the Osiris facility, with a 5 min period, lead¬ 
ing to an average thermal power of 66.5 MW. From the 
enthalpy balance performed online on the primary circuit 
of reactor we estimate a 2 % relative uncertainty HZI. 
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Furthermore the layout and burn-up of each assembly 
is made available, at the beginning of each reactor cycle. 
Detailed core simulations are then performed with the 
3D Monte-Carlo code TRIPOLI-4@ [TH]. We studied 
two core configurations, corresponding to the two typical 
initial burn-up maps of Osiris with either the control rod 
number 3 or 4 used as the last rod to control the reactor 
power. For each configuration two extreme positions of 
the control rods, at the beginning and end of the cycle, 
have been considered. During the cycle, the translation 
of the barycentre of the fissions due to the control rods is 
only 2.5 cm, so the deviation from the average barycentre 
of the fissions is negligible. The corresponding impact on 
the interaction rate is inferior to 0.3%. 

Between the two initial control rod sequences, the 
mean barycentre of the fissions is shifted by 3.6 cm at 
the beginning of a cycle and 1.1cm at the end, leading 
to a difference in the predicted rate inferior to 1 %. Still 
this effect is taken into account for the prediction of the 
mean neutrino rate. The overall uncertainty on the fis¬ 
sion fraction is assumed to be 2 %, covering the amplitude 
of the evolution during a full reactor cycle. Inserted in 
Eqjl] this leads to a 1 % uncertainty in the predicted neu¬ 
trino flux. The mean fission fraction are a235u = 92.6%, 
Q!239pu = 6.1%, q;238u = 0.8% uud Q!24ipu = 0.5%, as 
expected for a highly enriched fuel. 

The mean energies per fission Ek and the IBD cross- 
section are taken from m and m respectively. For 
the antineutrino fission spectra of each isotope we used 
the recently improved spectra from [2, converted from 
the ILL reference beta spectra of ^^®Pu and 

^■^^Pu [inilll]. The 238 u spectrum is taken from [i]. Con¬ 
sidering the dominant contribution from global 

2.2 % relative uncertainty is obtained for the mean inter¬ 
action rate of antineutrinos. 

The distance between the reactor core centre to the 
detector centre was determined using a survey of the 
geodesic combined with information from technical draw¬ 
ings of the reactor building. This led to a baseline of 
7.21 ± 0.11m, corresponding to a 3.1% uncertainty on 
the neutrino rate. The detector mechanical survey al¬ 
lowed us to compute the detector volume at 846.8 ± 7.0 L. 
Associated to the measurement of the target liquid mass 
of 739 ± 1 kg and the calculation of the ratio of free hy¬ 
drogen over carbon from the known chemical composition 
of the liquid (H/C=1.50 ± 0.02), we computed a proton 
density of 5.92 ± 0.05 • 10^® proton/m®. 

To propagate the Pg flux from the core to the Pg de¬ 
tected rate we developed a code called NuMC, with main 
task to integrate numerically Eq. and with a Monte- 
Carlo method. The code firstly computes the interac¬ 
tion rate for a perfect detector: a given fissile isotope is 
drawn according to the simulated fission fraction a^, then 
a Pg energy is drawn in the corresponding spectrum S'^, 
and finally a position in the core is drawn according to 
the simulated distribution of fissions in the core 'I'f,/c(^). 
From this stage, the Pg is propagated in a random direc¬ 
tion, and the length of the Pg path inside the detector 


(if any) is stored. Repeating a great number of times 
this procedure gives the average length of the Pg path 
inside the detector, which is Zgff weighted by the fission 
fraction. After numerical integration of the product of 
the spectra Sk and the cross-section cr^^, the interaction 
rate is then obtained by applying the correct normaliza¬ 
tion factor of Eq. We obtain for Osiris at 66.5 MW 
and 100 % of efficiency 910.8 Pg/day and 913.8 Pg/day for 
the two core configurations described previously, with a 
relative uncertainty of 4.6% (see table|l]). 

Then we can compute the expected number of Pg/day 
according to: 

= (i?™d3 X -h X X Cdet (3) 


with = 910.8 Pg/day (i?™'^^ = 913.8 Pg/day) the 

expected number of Pg interactions per day at 66.5 MW 
with the 3’'*^ (4***) control rod used as the last control fuel 
element during the cycle. = 0.32 and = 0.68 

are the relative amounts of lifetime in each configuration 
and Edet the detection efficiency (see section V F 11. We 
obtain = 913 Pg/day for Osiris at 66.5 MW and 

100 % of efficiency. 


Source 

Relative 
uncertainty (%) 

Baseline 

3.1 

Fission cross-section 

2.2 

Thermal power 

2.0 

Fission fractions 

1.0 

Number of protons 

1.2 

Total 

4.6 


TABLE I. Summary of relative uncertainties associated to 
the predicted Pe rate, with 100 % detection efficiency. 


In a second step, the NuMC code can be used as a 
Geant4 IBD event generator. An isotope and a fission 
vertex in the core are drawn, together with a vertex in 
the detector. This couple of vertices is stored or not ac¬ 
cording to the 1/L^ flux law. If this event is valid, an 
IBD event is generated according to the double differ¬ 
ential IBD cross-section m- Finally the full kinematic 
information of both the positron and the neutron are 
saved in a data file as an input for Geant4. 

We then simulate the Pg interactions in the detector, 
including scintillation and light collection on PMTs. The 
Geant4 output file is set to the same format as the Nucifer 
data. Finally processing our neutrino search algorithms 
on the simulated data determine the mean detection ef¬ 
ficiency, found to be 30.3% (see section VFl). 


IV. DETECTOR CALIBRATION 

The calibration system was designed to provide an ab¬ 
solute energy scale to the light response of the detector 
and to assess the linearity and the stability of the whole 
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detection system. Controlled light injections allow a cal¬ 
ibration from charge unit to photo-electrons (PE here¬ 
after). Radioactive sources deployed temporarily in the 
central tube inside the detector let us to perform the final 
calibration from PE-scale to the MeV-scale. 

The detector target vessel is equipped with 7 Teflon 
light diffusers, linked by optical hbres to Light-Emitting 
Diodes (LED) located in the electronic rack and con¬ 
trolled through the acquisition software. One of this LED 
is used to generate Single Photo-Electrons (SPEs) on the 
PMTs, 4 are used with different intensities to test the lin¬ 
earity and stability of the PMTs over their full dynamic 
range, and the last 2 could be used as spares. Sequences 
of 14 patterns, including 2 SPEs, random and multiple 
LED patterns, are continuously generated by the acqui¬ 
sition at a frequency of 5 Hz. What is referred to as a 
random pattern corresponds to no LED and, although 
periodic, is random with respect to physical triggers and 
therefore samples the pedestals of the digitization chain. 
This system allows the continuous measurement of both 
the pedestal and gain of each channel for each run, thanks 
to an automatic fit of the pedestals and of the SPE sig¬ 
nals. Each run and channel is consequently automatically 
calibrated in photo-electrons, allowing a proper sum of 
the 16 PMT charges. 



Total Charge (pe) 


FIG. 4. Detector response to a ®°Co source deployed at the 
centre of the target vessel compared to the detector response 
simulation. No threshold was applied in the simulation, lead¬ 
ing to the low energy discrepancy. For all calibrations, the 
hardware threshold is different from the data taking thresh¬ 
old, adapted to the energy range and the trigger rate. 

The detector response to multiple LED patterns com¬ 
pared to single LED pattern is linear within 1 %, show¬ 
ing no drift at higher intensities and the uncertainty on 
pedestals is below 1 %. On a run-to-run basis the vari¬ 
ance of the fitted SPE is typically 2.5% with perfect 
Gaussian shape. On longer time scales the gains and the 
linearity of the response to LED patterns have shown 
remarkable stability. Larger drifts due to temperature 
changes in the electronic rack are observed on the abso¬ 


lute value of pedestals and LED signals. They are corre¬ 
lated on all channels and corrected for each run. 

Small encapsulated radioactive sources can be de¬ 
ployed inside a vertical tube along the target central axis. 
Three gamma-emitting sources of few kBq activity have 
been used: a ^^^Cs source with one 661 keV y-ray, a ®°Co 
source with two y-rays of 1173 and 1332 keV in coinci¬ 
dence, and a ^^Na source with one y-ray of 1274 keV 
in coincidence with two y-rays of 511keV coming from 
the annihilation of the positron. We also used a neutron 
source of ^“^^Am-Be, emitting a neutron in coincidence 
with a y-ray of 4.4 MeV for ~ 75 % of the events. The 
activity of the ^^^Am reaches few MBq, leading to only 
~ 30 Hz of neutrons. This source is mainly used to test 
our analysis procedure for searching correlated pairs of 
events and study the neutron physics (capture time and 
detection efficiency). 

These calibration sources were inserted at different ele¬ 
vation levels in the central vertical tube and used to tune 
a Geant4.9.4-based simulation (in PE units) to reproduce 
the measurements. A good agreement between experi¬ 
mental and simulated spectra could be achieved for most 
positions (see Fig. |^. However at the highest source lo¬ 
cation, close to the acrylic buffer coupling to the PMTs, 
sizeable deviations from the data still remain, likely due 
to the difficulty of properly simulating light collection for 
events interacting close to the PMTs. 



FIG. 5. Detector response to a AmBe source deployed at 
the centre of the target vessel at three different periods from 
April 2014 to May 2015. The right peak corresponds to the 
neutron capture on Gd at 8 MeV. The middle peak comes 
from the reaction of a particle on ®Be, leading to the emission 
of a 4.4 MeV gamma ray and a fast neutron, and depositing 
a mean energy equivalent to an electron of 5.5 MeV. Due to 
the quenching effect, neutron induced nucleus recoils do not 
produce as much light as an electron of the same energy. The 
left peak corresponds to neutron capture on hydrogen, with 
the emission of a unique 2.2 MeV gamma ray. 

A global calibration factor of ~ 340PE/MeV has been 
obtained with the 4 sources of calibration. The uncer- 
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tainty on this factor is latter included in the global en¬ 
ergy scale uncertainty (see section VF2). With a source 
at the centre of the target we measure an intrinsic en¬ 
ergy resolution of 10 % at 1 MeV. This value increases to 
about 20 % for vertices uniformly distributed in the tar¬ 
get volume because of the light absorption on the Teflon- 
coated walls of detector target. Figure illustrates the 
very good stability of the detector response to the AmBe 
source along the whole data taking period. 


V. DATA ANALYSIS 

A. Data sample 


The signature of IBD events consists in a delayed coin¬ 
cidence starting by a prompt positron energy deposition, 
Epromptj followed by a neutron induced energy deposi¬ 
tion, Edeiay, due to the de-excitation gamma ray(s) after 
its capture on H or Gd within Ate+„. The principle of 
the analysis is to compare the number of detected 
to the prediction based on the reactor data and reactor 
core simulations. The challenge consists in the statistical 
separation of IBD events with respect to backgrounds in¬ 
duced either by random coincidences of single events or 
by correlated coincidences originating from air showers 
and possibly from reactor core radiations. 

At Osiris, a running cycle is nominally operating for 
three weeks with the reactor ON followed by a week of 
reactor OFF period. The results described in this paper 
are based on 10 cycles accumulated from June 2014 to 
July 2015, after the Nucifer upgrades. Power transients 
at the beginning of each cycle (few hours) are discarded 
from the data sample due to the difficulty to assess the 
thermal power of the reactor during this phase. There¬ 
fore our sample consists of 200 days of reactor ON at full 
power. Nevertheless the lifetime of the Nucifer data ac¬ 
quisition is only 145 days. A loss of 5 days is due to the 
100 ps veto of the acquisition after each tagged muons. 
The remaining inefficiency is dominated by unattended 
periods of data taking during which the synchronization 
of the different electronics readout buffers of the data ac¬ 
quisition system was lost. The bad runs are rejected by 
an automatic off-line quality check. The total of reactor 
OFF time is 153 days including inter-cycles periods and 
longer shutdowns for reactor maintenance. The Nucifer 
lifetime during this OFF period is 106 days. We thus ob¬ 
tain an overall data taking efficiency of 70 %. It is worth 
noting that this efficiency could easily be improved to 
more than 90 % by implementing an automatic recovery 
procedure to handle data acquisition system failures. 

In the following, event rates are reported with statisti¬ 
cal uncertainties only. Systematic uncertainties are gath¬ 
ered in the predicted neutrino rate uncertainty. 


B. Neutrino candidate selection 

The analysis cuts are optimized to reach the smallest 
relative uncertainty on the detected neutrino rate [22] in 
the severe background conditions described in the sec¬ 
tions IV Cl and IV D 11 

From the prediction of the emitted neutrino spectra 
and the kinematics of the IBD reaction the expected 
range of Eprompt is 1.022 MeV to about 7 MeV. Given 
the steep increase of background at low energy, a 2 MeV 
threshold is used, well above the hardware threshold, giv¬ 
ing negligible trigger inefficiency. Then the maximum 
prompt energy is set at 7.1 MeV above which the remain¬ 
ing IBD positron induced rate is negligible. 

The delayed event energy is not related to the incident 
neutrino energy. Its associated 8 MeV y-cascade from the 
neutron capture on a Gd nucleus is used to efficiently dis¬ 
criminate the neutrino signal against lower-energy back¬ 
ground events. However in a meter-scale detector like 
Nucifer one has to accommodate for important energy 
leakages shifting most of the delayed events into lower 
energy bins. Indeed, the most probable energy deposited 
in Nucifer after a neutron capture on Gd is about 4.5 MeV 
only. Again the limitation to open the Edeiay range to¬ 
wards lower energies is the steep increase of the acciden¬ 
tal background. We thus optimized the delayed energy 
range to 4.2 MeV < Edeiay < 9.6 MeV. 

The coincidence time window Ate+n is limited to 40 ps 
after each prompt candidate. This duration corresponds 
to only twice the expected neutron capture mean param¬ 
eter, T„. This upper limit is a compromise between the 
rejection of about 15 % of the neutrino signal and the 
mitigation of the accidental background rate, growing 
linearly with the prompt-delayed gate duration. Further¬ 
more the delayed event time must also be separated from 
the prompt energy deposition by more than 6ps due to 
the dead time of the QDC needed to process the PMT 
signals after each trigger. 


6 ps 

>fin IIS “id MS >60 US 



FIG. 6. Time selection of neutrino candidates. The prompt 
event must be separated by more than 100 ps with respect to 
any muon veto trigger. 

All events occurring less than 100 ps after a muon event 
are overridden. A muon event is defined either as a muon 
veto trigger or as the saturation of at least 15 PMTs. 
Finally an isolation selection criteria, called multiplicity 
cut, is applied to all prompt-delayed pairs. It imposes 
that no energy deposition occurs either 60 ps before the 
prompt event or 60 ps after the end of the delayed gate 














(see Fig. [^. This cut rejects cosmic-ray-induced back¬ 
ground with more than two particles in the same shower. 
All selection criteria are summarized in table El 


Griterion 

Applied cut 

Prompt event 

2.0 MeV < Eprompt < 7.1 MeV 

Delayed event 

4.2 MeV < Edeiay < 9.6 MeV 

Time selection 

6ps < Atg+„ < 40 ps 

Multiplicity 

No trigger 60 ps before prompt event 
No trigger 60 ps after delayed gate 

Muon veto 

100 ps 


TABLE II. Summary of neutrino selection criteria. 


In the following we will focus on the accidental and the 
correlated backgrounds. Indeed the very short baseline 
coupled with the shallow depth of the experiment imply 
that these backgrounds make much larger contributions 
than the expected neutrino signal does. Therefore an 
accurate determination of these spurious event rates is 
the main challenge of the Nucifer analysis. 


C. Accidental background 


When the reactor is OFF, the trigger rate is domi¬ 
nated by muons events 350Hz in the muon veto), 
and once muons have been excluded, by low energy nat¬ 
ural radioactivity decays (65.7 Hz above 2MeV) whose 
energy spectrum is shown in red on figure Beyond the 
rapid decrease of the energy spectrum after the thresh¬ 
old at about 500 pe, one can see a small peak at 2.6 MeV 
(880 pe) corresponding to the highest natural gamma ray 
of 208T1 decay (™Co in the tank stainless steel is less 
than 12 mBq/kg from a sample measurement and there¬ 
fore negligible), and an ankle near 8MeV (2700pe) at¬ 
tributed to natural neutron capture. The tail at high 
energy is attributed to the contamination by events in¬ 
duced by atmospheric showers in the vicinity to the de¬ 
tector that do not trigger the muon veto nor saturate at 
least 15 PMTs. 

When the reactor is ON the shape is different and the 
total event rate increase (177.1 Hz above 2 MeV). The as¬ 
sociated spectrum (blue curve on figure]^ shows a rather 
smooth shape up to about 10 MeV. Such high-energy de¬ 
positions are attributed to neutron radiative captures on 
metallic structures and concrete wall components, such 
as iron, nickel or aluminium, occurring in the vicinity of 
the detector and emitting high-energy y that can pass 
through the shielding and finally reach the detector tar¬ 
get. The bump around 6 MeV (2000 pe) comes from the 
decay gamma rays in the deactivation circuit (see 
section 


HB). 


As a consequence both the prompt and the delayed en¬ 
ergy windows are affected by the reactor induced events, 
which has a strong consequence: the accidental back¬ 
ground event rate scales quadratically (instead of lin¬ 
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FIG. 7. Energy spectrum of the single event rate at reactor 
full power and during reactor OFF combined periods. 


early) with the gamma-ray flux 4>y^ proportional to the 
reactor power. A crude model of the neutrino-like event 
rate is Tp = Tprompt x Tjeiay X Atg+„. Both prompt and de¬ 
layed event rates are dominated by backgrounds events, 
and if both are affected by reactor induced events we get 
Tprompt OC (j)y , Tdelay oc (j)y and finally Tp (X Adding 
lead shielding could therefore reduce greatly the acciden¬ 
tal background, but the mass limit on the Osiris floor has 
already been reached. 



FIG. 8. Distribution of the time difference between prompt 
and delayed events for all correlated candidates (blue) and for 
the accidentals (red). The data sample accumulates 10 reac¬ 
tor cycles. 

Nevertheless this high accidental background level can 
be accurately measured by analysing off-time coinci¬ 
dences. For each delayed-like event, all past prompt 
events in the run are shifted in time up to 100 times by 
steps of 1 ms, until the next past prompt is repositioned 
less than 1 ms before the current delayed event. If the so- 
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formed virtual pair complies with the neutrino selection 
cuts, it is counted as an accidental event. This method 
allows to measure the accidental background 100 times 
for each run, in the exact same conditions of data tak¬ 
ing, and consequently this background is measured with 
a negligible statistical uncertainty. Unfortunately, this 
method limits the impact of the high accidental event 
rate on the final statistical uncertainty only for one of 
the two terms of the subtraction (see also table [V|). Not¬ 
ing N the event number. At the acquisition time, and r 
the events rates, the correlated event rate error is: 



where we see that if the 100 shifted gates allow to ne¬ 
glect the uncertainty brought by the accidental events 
subtraction, the final uncertainty is still dominated by 
the accidental event rate, as Tcorr ^ Tacc- 

Furthermore, great care is brought into the evaluation 
of the efficiencies of the multiplicity cut and muon veto 
as they slightly differ between the searches for correlated 
pairs and accidental pairs. A cross-check of our accurate 
determination of accidental rate is illustrated in figure 
where the number of measured correlated and accidental 
pairs is plotted as a function of Ate+n for a same sam¬ 
ple of reactor ON data. The correlated events clearly 
show up as an exponential curve on top of the acciden¬ 
tal background with a decay time corresponding to the 
n-capture time r„ ~ 20 ps. The distribution of pure ac¬ 
cidentals, determined with the off-time prompt events, 
perfectly matches the correlated curve for Atg+„ t„. 

Within the selection cuts, summarized in table [IT] the 
average accidental rates measured for the whole data 
sample are 69.1 ± 0.1 events/day when the reactor is 
OFF and 3476.3 ± 0.7events/day when the reactor is ON 
(at full power), respectively. Therefore this dominating 
background is 11.9 times larger than the expected neu¬ 
trino rate. In the following all quoted rates will be cor¬ 
rected for their associated accidental rates. 


D. Cosmic induced background 

1. Origin and rejection 

The correlated events above the accidentals events in 
figure 1^ are still not a pure sample of neutrinos. At the 
shallow depth of Nucifer (12mwe), the cosmic-ray parti¬ 
cles induce correlated backgrounds through the multiple 
secondary particles produced in the air shower and also 
through their interaction in the ceilings above the detec¬ 
tor (especially for the penetrating muons). It is worth 
noting that the overburden above the detector is not 
sufficient to stop all the hadronic component of the at¬ 


mospheric cascades, leaving the possibility for some fast 
neutrons to reach the liquid scintillator. 

The most probable candidates at the origin of cosmic 
induced correlated backgrounds are fast neutrons from 
a air shower or created by an inelastic muon interaction 
(spallation) in materials above or nearby the detector. 
First, a fast neutron can scatter off nuclei in the detec¬ 
tor target, mimicking prompt-like energy deposition, and 
later be captured on Gd providing a delayed-like energy 
deposition with a similar prompt-delayed time correla¬ 
tion than expected for neutrinos. Second, two neutrons 
from the same shower can be captured successively in 
the detector, after some diffusion in the liquid or in the 
shielding. As the gamma collection efficiency after neu¬ 
tron capture on Gd is poor in Nucifer, the first neutron 
capture energy deposit is very likely to be in the prompt 
energy window. Third, fast neutron can also produce 
high energy gamma rays by inelastic scattering on nu¬ 
clei, such as the ^^C{n,n')y reaction producing 4.4 MeV 
gamma rays on the first nuclear excitation level. The 
gamma interaction would of course mimic the prompt 
event and the later neutron capture the delayed event. 

Fortunately, this last background process whose signa¬ 
ture is very similar to the IBD can be safely neglected 
with respect to the first process (elastic scattering mim¬ 
icking the prompt event): the ^^G(n,n')Y cross-section 
is always well below both the hydrogen and the carbon 
elastic scattering cross-section, by a factor 5 to 10 each 
depending on the energy El- 

Highly energetic muons may also create long-lived |3-n 
emitters such as ®Li or ®He when interacting with Garbon 
nuclei belonging to the scintillator molecules. The rate 
of (3-n events from ®Li or ®He decays was estimated to 
2.7 events/day by scaling the measurement of the Dou¬ 
ble Ghooz experiment to the Nucifer shallow depth |23j . 
We compared this estimation with our analysis of events 
following high energetic muon depositions in the detec¬ 
tor target at time differences compatible with the decay 
time of ®Li (~ 270 ms). The muons showering inside the 
Nucifer target could be tagged by lowering the gain of 
one target PMT. No ®Li induced background candidate 
was clearly identified. In consequence an upper limit of 
less than 12events/day at 95% C.L. was set [52]. There¬ 
fore this background can be considered as negligible with 
respect to the total rate of correlated events. 

As the muons entering the detector are tagged by the 
veto system surrounding the detector, most of the sec¬ 
ondary products of the muon interactions are removed 
from the data sample by discarding events occurring 
within 0.1ms after each recorded muon. A further back¬ 
ground reduction is provided by the multiplicity cut, re¬ 
jecting most of the neutron induced background. In fine, 
the remaining contribution is subtracted using reactor- 
OFF data, representing about 45% of the data taking 
time. Averaged over the whole reactor OFF data, a corre¬ 
lated event rate of 1145.4 ± 3.4 ± 2.5events/day is mea¬ 
sured, i.e. 3.9 times higher than the expected neutrino 
rate. 
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2. Stability and correction 

Last but not least, before subtracting this correlated 
background to obtain the experimental neutrino rate, one 
has to assess its stability over the time scale of a reac¬ 
tor cycle. Indeed the correlated background is found to 
vary over time and consequently an extra probe of the 
correlated background must be defined to monitor and 
subtract properly the actual background rate during each 
reactor ON period. 

A probe is the muon rate, accurately recorded for each 
run. Here we define a muon trigger as the combination of 
a signal in the muon veto and the saturation of at least 
15 PMTs. This definition reduces the sensitivity to high- 
energy y that can trigger the muon veto when the reactor 
is ON and allows a reliable comparison between reactor 
ON/OFF muon rates. Table [V| shows that indeed the 
muon rates measured this way for the averaged reactor 
ON and OFF periods are very close, 110.6 and 110.1 Hz 
respectively. The main origin of this variation of muon 
rates is the natural variation of the atmospheric pres¬ 
sure: we could check that the remaining 0.5 Hz difference 
is coherent with the different mean atmospheric pressures 
recorded for the two averaged periods (998 and 1002 hPa) 
and the observed correlation of —0.117 ± 0.008 Hz/hPa. 
Figure]^ (top) shows that the variations of the correlated 
event rate among various reactor OFF periods are very 
well correlated with the muon rate. Once the correla¬ 
tion with pressure is removed no sizeable dependence on 
temperature or humidity was found. 




FIG. 9. Correlation between the rate of corre¬ 

lated background measured with the reactor ON/OFF and 
the rate of triggers tagged as muons. The slope is 
—22.8 ± 5.0 event day“^ Hz“^. 

At first glance one would expect a zero intercept and 
a relative change of the muon rate equal to the relative 
change of the correlated background. This is clearly not 
observed. More complex physics must be brought into 
play to explain the true dependence, probably due to the 


fact that the tagged muons saturating the detector are 
not fully representative of those inducing the correlated 
background (coming from interactions outside the detec¬ 
tor) . Hence a linear fit is considered more as an effective 
model within the limited range of variation of the muon 
rate. What matters for the subtraction is that the same 
correlation applies on the reactor-OFF data. This is in¬ 
deed the case as shown in figure (bottom) where the 
muon rate dependence measured for reactor OFF periods 
is in very good agreement with the one for reactor ON 
data. Hence the subtraction of the correlated background 
is performed the following way: 

- The absolute correlated event rate per day is given 
by the reactor OFF data. The global correlation 
with the muon rate is determined as accurately as 
possible, combining reactor ON and OFF data to 
constrain the muon rate dependence. The final er¬ 
ror band is shown as a grey shaded area in figure 

- The muon rate is measured for each reactor ON 
period. 

- The corresponding correlated background and its 
associated uncertainty is deduced from the corre¬ 
lation with the muon rate. This reactor OFF con¬ 
tribution is subtracted from the total rate of corre¬ 
lated events measured reactor ON. 

As the difference between average reactor ON and OFF 
muon rate is only 0.5 Hz and the slope of the correction 
is —22.8 ± 5.0 event day“^ Hz“^, the correction reaches 
11.4 event/day and introduces a small systematic uncer¬ 
tainty on the final neutrino rate of 2.5 event/day. When 
the data sample is separated in shorter time period, this 
error increases for each single data point, depending on 
the actual muon rate. 


E. Reactor induced correlated background 


1. Neutron induced correlated background 


Due to the proximity of the detector with the reac¬ 
tor, fast neutrons generated by the core could penetrate 
the polyethylene shielding producing additional corre¬ 
lated background. This kind of events could be even 
more problematic with respect to the cosmic ray induced 
neutrons since the associated background cannot be sub¬ 
tracted using reactor OFF data. We consider here only 
one correlated background process: neutron elastic scat¬ 
tering on nuclei mimicking the prompt event. Even if 
a fission generates several fast neutrons, it is extremely 
unlikely that 2 neutrons from the same fission scatter 
off nuclei in the (small) Nucifer detector located several 
meters away. We also consider that the inelastic scatter¬ 
ing on nuclei process is negligi ble wit h respect to elastic 
neutron scattering (see section VD 1). 


In order to check if a background contribution remains 
on top on the reactor OFF data, we exploit the Pulse 
Shape Discrimination (PSD) capability of the Nucifer liq¬ 
uid scintillator. Aside the total charge recorded for each 
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trigger (Qtot), a late charge is also recorded (Qtaii, see[n|). 
The distribution of the ratio Qtaii/Qtot of all prompt can¬ 
didates is displayed in figure [l0| for reactor ON and OFF 
data. As expected, two peaks emerge, one centred at 
Qtaii/Qtot — 0.24 corresponding to electron recoil and an¬ 
other one centred at a higher value, 0.29, corresponding 
to nucleus recoils induced by fast neutrons. In the com¬ 
parison between the ON and OFF PSD spectra the neu¬ 
trino signal clearly shows up as a Gaussian distribution 
centred at 0.24, as expected from the positron (electron¬ 
like) recoil. For large Qtaii/Qtot values corresponding 
to nucleus recoils, both spectra overlap and demonstrate 
that the background subtraction is well under control. 

Adjusting each curve with the sum of two Gaussians, 
the compatibility of the nucleus recoil peaks in our errors 
bars is confirmed. However, the reactor ON peak ex¬ 
hibits 18 more events than the reactor OFF peak, in the 
0.3 to 0.4 Qtaii/Qtot range. We conclude that the reactor 
induced correlated event rate is 0 ± 18 events/day, and 
add this as an additional systematic uncertainty to the 
final neutrino rate measurement. It is worth noting that 
the delayed gate used to record Qtaii has suffered from an 
electronic jitter during the major part of our data taking, 
not constant in time, before we could identify and cor¬ 
rect the defective gate generator. This electronic jitter 
is our best hypothesis to explain the mismatch between 
the nuclei recoil peaks ON and OFF. 

To better constrain the neutron induced correlated 
background, we performed a complete TRIPOLI-4@ 
simulation. The starting point is a core simulation that 
samples outgoing fast neutrons (> 2MeV) at the surface 
of the core. These particles are then used as a source 
for a second simulation that propagates the neutrons to 
the Nucifer target volume, using sophisticated variance 
reduction techniques to push the particles through the 
thick absorbing media. In particular, we used INIPOND, 
a special built-in module of TRIPOLI-4@ based on the 
exponential transform method |24j . This exponential bi¬ 
asing is performed using an importance map which pro¬ 
vides information on the probability, for each point of 
the phase space, for a particle to reach the detector. It 
is first calculated with a simplified deterministic solver 
and then adjusted by hand. The code uses this informa¬ 
tion to adjust the propagation of the neutrons along the 
path to the detector, thereby reducing the variance of the 
calculated observables without introducing any bias. 

We found that the neutron elastic scattering rate 
on hydrogen in the whole Nucifer target volume is 
4 • 10“® event/day for neutrons with energies higher than 
2 MeV, corresponding to an attenuation of the order of 
10“^^. This rate is a large overestimation because we 
neglected the real energy transfer, the quenching and 
the prompt-delayed coincidence, considering each neu¬ 
tron elastic scattering on hydrogen as a background event 
satisfying all our selection cuts. The statistical uncer¬ 
tainty reported by TRIPOLI-4@ is at the percent level 
but large systematic uncertainties have to be considered, 
simply by the propagation of errors on initial parame¬ 


ters (cross-section, composition...) on such enormous 
attenuation. The results should therefore be considered 
as an order of magnitude estimate. Despite the large 
uncertainty, our conservative approach ensures that the 
neutron induced correlated background is negligible. 



FIG. 10. Distributions of the PSD parameter Qtaii/Qtot for 
reactor OFF (squares) and reactor ON data (circles), obtained 
after subtraction of the accidental background and correction 
for the cosmogenic rate modulation. 

Note that the separation between proton and electron¬ 
like recoils could in principle be used for further rejection 
of the correlated background. A cut Qtaii/Qtot < 0.27 
would remove about 2/3 of the reactor OFF background. 
However this criterion is not currently in force in our 
analysis since the rather poor separation between the 
2 peaks implies a potentially large migration of events 
from both sides of the cut should slight drifts in the pulse 
integration occur. The application of this PSD cut would 
thus lead to an extra systematic uncertainty whereas our 
dominant background is coming from accidental coinci¬ 
dences. 


2. Gamma induced correlated background 

We also considered the gamma correlated background: 
a gamma ray with sufficiently high energy could excite 
a nucleus by inelastic scattering, and if excited above 
the neutron separation energy the nucleus could emit a 
neutron. This process demands very energetic gamma 
ray since the typical neutron separation energy is above 
10 MeV for stable light nuclei (^^G, ^®0, ®®Fe) and in the 
6-8 MeV range for heavy nuclei (Pb, Gd), with strong 
effects of neutron number parity (several MeV) so that 
some rare stable nuclei have lower neutron energy sepa¬ 
ration (e.g. 4.1 MeV for 4.9MeV for ^^C). 

If the gamma energy is high enough, the neutron could 
be fast and create correlated background as in the previ¬ 
ous scenario. But up to 6 MeV proton recoils produce 
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3 times less light than electron recoil due to quench¬ 
ing [25l [26], therefore the neutron should have at least 

6 MeV to enter our prompt energy window. Adding this 
energy to the neutron separation energy demands an ini¬ 
tial gamma energy higher than what is observed in the 
Nucifer spectrum (see Fig. [^. 

Another possibility is the creation of a low energy neu¬ 
tron, the remaining energy being left to the gamma ray. 
If both particles enter the detector and are absorbed, the 
gamma ray can create the prompt event and the neu¬ 
tron capture the delayed event, with a time correlation 
identical to neutrino events. To enter the prompt energy 
windows, the gamma ray must have kept at least 2 MeV, 
and because the neutron would have low energy, the prob¬ 
ability to go through the polyethylene shielding without 
being captured on Boron would be negligible. There¬ 
fore the initial inelastic scattering should take place in 
the tank, on some odd neutron number light nuclei. But 
such nuclei have low natural concentration (percent or 
less), at such energy the cross-section of (y,n) reaction 
(2.0 • 10“^ barn for at 7 MeV [13 [55]) is much lower 
than Compton scattering cross-section (3.96 • 10“^ barn 
for C at 7MeV [29]) and the interaction rate of high en¬ 
ergy gamma in Nucifer is low enough (about 0.5 Hz above 

7 MeV) so that this process can be safely neglected. 


F. Detection efficiency and associated systematics 

1. Detection efficiency 

As explained in section]!!] a dedicated software package 
generates a set of simulated Vf, events in the same format 
than the data. To determine the detection efficiency we 
apply the exact same analysis chain and compute the 
amount of rejected events. Since the path lengths of y- 
rays and neutrons in the liquid scintillator (^ 20 cm) are 
not negligible with respect to the size of the detector 
target (1.2m diameter), edge effects induce sizeable cor¬ 
relations between the analysis selection criteria. In con¬ 
sequence the total detection efficiency is not the simple 
product of the efficiencies of all single cuts. It is defined 
as the ratio of the number of simulated events passing all 
the analysis selection criteria to the total number of neu¬ 
trino vertices generated in the liquid scintillator. Using 
the selection procedure described previously the global 
detection efficiency is found to be 30.3%. 

As shown on table [III] taking each cut individually the 
selection of the delayed energy deposition (47.9% effi¬ 
ciency) is the selection criteria that degrades most the 
efficiency. This is due to the lack of containment of 
the multiple high-energy y-rays released after the neu¬ 
tron capture on Gd in the small, 619.2 mm in radius 
and 704.9 mm in height, detector target. The time selec¬ 
tion is the next limiting cut (63.1 %), rejecting the same 
amount of event below and after the gate, and finally the 
low energy cut at 2 MeV has a rather strong influence 
on the prompt event selection (76.6%). According to 


the simulation, IBD neutrons are captured mostly by the 
Gadolinium (85%), then by the Hydrogen (11%), the 
Boron (2%, in the polyethylene shielding) and the re¬ 
maining by the steel components. 7.2 % of IBD neutrons 
are captured outside the target. 

The efficiency of the multiplicity cut is computed by 
counting the number of diode triggers with and without 
applying this cut. Because of their different single rates 
reactor ON and OFF have different multiplicity efficien¬ 
cies found to be 0.974 and 0.990, respectively. They both 
have a negligible associated systematic uncertainty and 
the correction is applied during the data analysis. Fi¬ 
nally the 3.9 % dead-time induced by the p.-veto is inde¬ 
pendently computed with negligible uncertainty from the 
p. triggers, taking in account overlaps between veto gates. 
This veto time is subtracted from the total live-time of 
the analysed runs. 


Selection cut 

Efficiency 

6 ps < At < 40 ps 

63.1% 

2.0 MeV < Eprompt 7.1 MeV 

76.6 % 

4.2 MeV < Edeiay < 9.6 MeV 

47.9 % 

Cdet 

30.3 % 

Dead-time 

Multiplicity ON 

2.6% 

Multiplicity OFF 

1.0% 

Muon veto 

3.9% 


TABLE III. Detection efficiency or induced dead-time of the 
selection cuts. As explained in the text the total efficiency is 
not the simple product of cut efficiencies because of correla¬ 
tions between the cuts. 


2. Systematic uncertainties 

The uncertainty on the energy scale is determined by 
comparing the delayed spectrum of correlated pairs in 
the reactor OFF data (Fig. ]TT] ) with the Geant4.9.4 sim¬ 
ulation of thermal neutron uniformly distributed inside 
the lead shielding, once the Geant4 simulation has been 
tuned to reproduce at best the source spectra at differ¬ 
ent positions. This set of data was chosen because it 
is independent from the calibration data, it shows sev¬ 
eral spectral features among the whole energy range and 
the energy depositions are uniformly distributed in the 
whole target volume, as expected from neutrino candi¬ 
dates. Note that the source in the simulation, uniformly 
distributed thermal neutron, is a crude approximation of 
the natural neutron flux at shallow depth. But the simu¬ 
lation is not intended to reproduce exactly the measured 
spectrum, particularly the relative weight of the energy 
features, but to allow the comparison of the position of 
these features to establish the energy scale uncertainties. 

The main features of the energy spectrum (n-capture 
peak on H at 750 PE, a middle bump at about 1700 PE 
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and a high energy edge of n-capture on Gd) are clearly 
visible for both data and simulation. The n(H,D)y 
peak position, the transition between the two plateaux 
at 1200 PE and 2200 PE, and the n(Gd, Gd)y shoulder 
position are fitted on both data and MG. The largest 
difference is found to be 2.5% and is taken as a safe es¬ 
timate of the energy scale uncertainty. To determine the 
efficiency uncertainty due to energy scale, we repeated 
our analysis of the simulated neutrino data set 1000 times 
with a variation of energy cuts. At the beginning of each 
analysis, a factor is shot in a Gaussian distribution of 
relative standard deviation 2.5 %, and each energy cut is 
multiplied by this factor. The distribution of all the neu¬ 
trino rates found is a Gaussian which standard relative 
deviation, 1.0%, is the uncertainty on efficiency due to 
energy scale. 

Nucifer is made of a unique volume of liquid scintillator 
enclosed in a stainless steel tank. Therefore, neutrino 
interaction outside the target leading to light production 
can only happen in the acrylic buffer, with the positron 
depositing at least 2 MeV in the liquid scintillator. Thus 
only the first centimetres of the buffer are concerned. 
We simulated IBD events uniformly distributed in the 
buffer and found that less than 0.15% of such events 
pass the analysis cuts, a negligible quantity (less than 
0.5event/day). 



FIG. 11. Photo-electron spectrum of delayed energy events 
during a reactor OFF run compared to a Geant4.9.4 simula¬ 
tion of low energy neutrons uniformly distributed in a sphere 
of 2 m radius centred on the target vessel. The n-capture peak 
at 750 PE (2.2 MeV) is prominent because of the numerous 
captures occurring in the polyethylene shielding around the 
target. The data threshold comes from the software analysis, 
thus not affecting the n(H, D)y peak, while the simulation has 
no threshold. 

The uncertainty associated to the time selection crite¬ 
ria is determined by comparing Geant4.9.4, Geant4.9.4 
with the addition of a module dedicated to thermal neu¬ 
tron physics [3D] and TRIPOLI-4 (r) simulations, using 
different description of low energy neutron physics. We 


find a 1.7% error from the dispersion of these 3 sim¬ 
ulation results, taking as central value the Geant4.9.4 
simulation with the additional thermal neutron module. 
We also add a 1 % uncertainty on the efficiency due to 
the modelling of the gamma cascade following a neutron 
capture on Gd, by comparing the results of our 3 simu¬ 
lations. 

The measurement of the neutron mean capture time 
with the AmBe source or the reactor OFF correlated 
events can not be directly used because the initial energy 
and spatial distributions of the neutrons are different. 
Indeed, the thermalization phase of neutrons from IBD 
is not detectable due to our electronic dead time (6ps), 
while the Aty-n distribution of AmBe events shows a 
peak at 10 ps, then an exponential decay with a shorter 
parameter of 17.6 ± 1.2 ps for a source at the centre of 
the detector. Neutrons from a^^^Cf source would suffer 
from the same problem. 

Finally, the efficiency and the associated uncertainty 
are 30.3 ± 2.2 % (see table IVI, but here 2.2 % is an ab¬ 
solute uncertainty on the efficiency of 30.3 %. Therefore, 
the relative uncertainty on efficiency is 2.2%/30.3% = 
7.2%, and the predicted neutrino detection rate is 
= 9I3Te/day x 30.3 ±2.2% = 277±2 3Pe/d ay 
(see section for the predicted interaction rate IIIB). 


Source 

Uncertainty 

Time cut 

1.7% 

Energy scale 

1.0% 

Gd(n,y) cascade 

1.0% 

IBD in buffer 

0.2% 

Absolute efficiency uncertainty 

2.2% 

Relative efficiency uncertainty 

7.2% 

z7e interaction rate 

4.6% 

Total prediction uncertainty 

8.5% 


TABLE IV. Breakdown of the detection efficiency systematic 
uncertainties, and final prediction uncertainty. Pe interaction 
rate uncertainties are listed on table [H 


VI. RESULTS 

The distribution of neutrino events is obtained after 
subtracting the accidental background and the correlated 
background from the raw number of candidate pairs. In 
Fig. I^the spectrum of the prompt-delayed time differ¬ 
ences is found to be in agreement with a pure exponential 
curve with a decay time of 19.7 ± 0.9 ps, compatible with 
the decay time expected from the simulation. 

We can use this curve to set a constraint on possi¬ 
ble residual accidental background that would essentially 
show up as an extra offset. Adding this free parameter in 
the fit the value converges on 8 • 10“^ with an uncertainty 
of 6.6 • 10“^ that makes it perfectly compatible with zero. 
Integrating this error bar between 6 and 40 ps we obtain a 
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FIG. 12. Distribution of the prompt-delayed time differences 
for all neutrino candidates after background subtraction. The 
spectrum of the prompt-delayed time differences is found to 
be compatible with a pure exponential curve with a decay 
time of 19.7 ± 0.9 ps, in agreement with the expected value. 

possible “accidental-like” contribution of 2.3events/day. 
This corresponds to only 0.07 % of the measured acciden¬ 
tal rate and thus demonstrates the good quality of the 
subtraction. 

From the number of entries we get a total of 40 760 Fe 
detected in Nucifer, corresponding to a mean rate of 
^obs _ 281 ± 7(stat) ± 18(syst) Fe/day. The expected 
detected neutrino rate is = 277±23Fe/day and 

therefore = 1.014 ± 0.108. 


evolution of the detected neutrino rate is shown in fig¬ 
ure where the data are grouped in periods of about 
5 days. The alternation of ON-OFF periods is clearly 
visible. The null rate corresponds to the correlated back¬ 
ground event rate averaged over the full reactor OFF 
data. 

The statistical accuracy per 5 days of data taking is 
about 11 %. As expected the intrinsic statistical accuracy 
of the daily neutrino rate is decreased by a factor 4 due to 
the total background contribution being 16 times higher 
than the signal (11.9 from the accidentals and 3.9 from 
the cosmic-ray-induced correlated events). As illustrated 
in figure the isotopic composition of the Osiris reac¬ 
tor is quite stable in nominal operation. Consequently 
only sub-percent variation of the emitted Fg flux are ex¬ 
pected between the beginning and the end of a cycle m- 
The observation of this Osiris intra-cycle rate evolution 
is clearly out of reach of a Nucifer-like detector. 

The current statistical sample accumulated by the Nu¬ 
cifer experiment provides a modest sensitivity to test the 
Reactor Antineutrino Anomaly. The current 7.2 % large 
uncertainty on the absolute predicted normalization fac¬ 
tor prevents any strong conclusion in favour or against 
the averaged discrepancy, as shown by fig¬ 

ure Improving these results is beyond the scope of 
this publication. A significant work will be necessary to 
refine the prediction of the expected neutrino rate, in¬ 
cluding an improved determination of efhciency, and to 
reduce the reactor induced background systematics. 


VII. SENSITIVITY TO THE PLUTONIUM 
CONTENT OF THE CORE 


Single rates (Hz) 

ON 

OFF 

|r veto 

344.3 

376.4 

q veto & Saturation 

110.6 

110.1 

Singles 

177.1 

65.7 

Prompt singles 

75.4 

16.1 

Delayed singles 

15.7 

1.6 

Pairs (/day) 

ON 

OFF 

Candidates 

4903 ± 7 

1223.5 ±3.4 

Accidentals 

3476.3 ± 0.7 

69.1 ±0.1 

Correlated 

1426±7± 18 

1145.4 ±3.4 ±2.5 

Dobs 

281 ± 7(stat) 

± 18(syst) Pe/day 


277 ± 

23Pe/day 


TABLE V. Summary of relevant rates for reactor ON and 
OFF periods. Note that the number of observed neutri¬ 
nos is not directly the ON-OFF difference of correlated 
events, as correction for the multiplicity cut and for the dif¬ 
ferent mean muon rates have been applied. Errors reported 
on candidate and accidental event rates are statistical only, 
while errors on correlated rates also show systematic uncer¬ 
tainties due to background subtraction. All other systematic 
uncertainties are applied on the predicted rate. 

All relevant rate, averaged over the full reactor ON 
and OFF data sets, are summarized in table [V] The time 


With the end of the Cold War, hundreds of tons of 
weapon-grade plutonium were determined to be surplus 
to U.S. and Russian defence needs. In April 2010 the US 
and Russian governments signed a protocol amending the 
2000 Plutonium Management and Disposition Agreement 
(PMDA), which commits each country to dispose of no 
less than 34 metric tons (MT) of excess weapon-grade 
plutonium and envisions disposition of more weapon- 
grade plutonium over time. The combined amount, 
68 MT, represents enough material for several thousands 
of nuclear weapons. The current approach is to trans¬ 
form the weapon-grade plutonium into mixed oxide fuel 
and irradiate it in reactors. We study here how a small 
neutrino detector like Nucifer could monitor this kind 
of operation to guarantee that plutonium is really being 
burnt in the reactor. 

The neutrino rate measured by Nucifer is used as a 
calibration point for a nuclear fuel highly enriched in 
235u with 92 % of the fissions coming from this isotope 
(see section III). This approach can be seen as a mea¬ 
surement of our absolute normalization with the statisti¬ 
cal 2.5 % accuracy (from 281 ± 7(stat) ± 18(syst) Pe/day) 
provided by our sample. 

Using the MCNP Utility for Reactor Evolution 
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Days 


FIG. 13. Antineutrino rate measurement monitoring the Osiris nuclear reactor operations. Each data point (blue diamonds) 
corresponds to about 5 days of data taking with its associated statistical uncertainty. The grey shaded area is the rate 
expectation above the mean correlated background when the reactor is OFF, referred as zero level here and plotted as red dots. 
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FIG. 14. Comparison of the Nucifer measurement of 
with others short baseline experiments and the 
average value of the RAA [5]. Without Nucifer, the average 
ratio is 0.938 ± 0.024 (2.6 cj), and with Nucifer the average is 
0.940 ± 0.024 (2.5 ff). 


(MURE) [32 [33] , we were able to simulate the evolution 
of the Osiris fuel with a full 3D Monte-Carlo simulation. 
Figure shows an example of the evolution of the fission 
fractions of ^^®Pu, and ^'^^Pu as predicted by 

MURE [3T] . Starting from a fresh uranium fuel enriched 
at 19.75% an equilibrium is reached after the 7*^ cycle. 
As expected the fission of the dominates by far and 
all fission fraction evolutions during a cycle are moderate. 

We then simulated the operation of the Osiris reac¬ 
tor with part of its fuel elements loaded with ^^®Pu, like 
in UOX-MOX fuel. The MURE simulation of the fuel 
evolution being very time consuming, we used for this 
purpose a simple model of the core [34j : we solved the 



FIG. 15. Fission fraction in % of each hssile isotopes, 

238-y^ 239p^^ as predicted by the MURE simulation 

(solid lines) jST]. Here the core is initially filled with fresh 
fuel and its isotopic composition evolves for 14 cycles, but 
all actual Osiris cycles correspond to the equilibrium regime. 
The dotted lines illustrate the effective approach dedicated to 
study the sensitivity to the ^^®Pu mass in the core after it is 
tuned to reproduce the reference MURE results. 


Bateman’s equations of the evolution of the uranium and 
plutonium isotopes assuming a uniform neutron flux ad¬ 
justed to match the mean Osiris thermal power and with 
effective nuclear cross sections (n-capture and fission) cal- 
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ibrated to reproduce the MURE simulation results. Fig¬ 
ure [T^ shows that a good agreement could be reached for 


the nominal core composition of Osiris, good enough for 
the purpose of this study. 



As expected the neutrino rate decreases as the mass 
of ^^®Pu in the core, hence the contribution of ^^®Pu fis¬ 
sions, increases. This is illustrated in figure where 
a global normalization factor is applied to our model in 
order to set the prediction for the nominal core compo¬ 
sition equal to the observed rate. Then assuming the 
same data taking period and background conditions, the 
detection of a lower Pe rate with 95 % confidence level 
is equivalent to a drop of 2.23 times of our current rate 
uncertainty (7.1), i.e. 15.8 less detected per day. This 
sensitivity limit is reached for about 1.5 kg of ^^®Pu in the 
core, a mass representing 10 % of the total mass of fission¬ 
ing elements in the Osiris core at equilibrium. As already 
stated before, this sensitivity is driven by the large acci¬ 
dental background at the Osiris site. If the same Nucifer 
detector was installed close to a commercial reactor with 
a baseline (25 m) and background conditions (S/B=4) 
similar to the SONGS experiment [35], the same sensi¬ 
tivity could be reached in less than 3 days. 


VIII. CONCLUSION AND OUTLOOK 


FIG. 16. Mean neutrino rate detected by Nucifer (blue 
diamonds), concatenating Osiris cycles. The grey area shows 
the interval at 95 % C.L. The other horizontal lines illustrate 
the predicted evolution of the detected flux for an increasing 
mass of Plutonium in the Osiris core. 

In standard operation the Osiris core is loaded with 
14.00 ± 0.75 kg of 235U and 450 ± 50g of ^a^Pu, the er¬ 
ror bar representing the variation of the masses around 
these mean values during a 21 days cycle. This equilib¬ 
rium is reached after an initial load of about 20 kg of 
235u_ rpp study the impact of plutonium content on the 
detected neutrino flux we tried various initial configura¬ 
tions with up to 10 kg of Pu in the initial composition 
of the core. However the initial mass of all fissile mate¬ 
rials (mainly 235^ ai^yj^ays kept equal to 

20 kg for all configurations. Then we simulated the reac¬ 
tor core evolution during several cycles at constant full 
power, until the equilibrium was reached. At this stage 
we determined the associated mean mass of in the 
core from the number of atoms predicted by our model. 
Finally, the mean neutrino rate was computed from the 
number of fissions per isotope in our model and from the 
reference ratios of detected per figgion listed in 
table rvn 


Isotope 

235u' 

23S-y 

239pg 241pg 

Relative Vf. rate 

1 

1.512 

0.635 0.900 


TABLE VI. Variation of detected antineutrino rate per fission 
normalized to one fission of The reference antineutrino 

spectra per fission are taken from [3] and [4] and interaction 
cross-section from m- The prompt energy cuts of the Nu¬ 
cifer analysis are applied corresponding to 2.78 MeV < < 

7.88 MeV. 


This article reported on the features and performances 
of the Nucifer experiment operating at the Saclay re¬ 
search centre of the French Alternative Energies and 
Atomic Energy Commission (CEA) since Spring 2012. 
The experimental configuration, the Osiris research nu¬ 
clear reactor, and the detector setup have been presented. 

We discussed the installation and the operation of the 
detector 7.21 ± 0.11 m away from the Osiris reactor, mak¬ 
ing Nucifer the second world-shortest baseline neutrino 
experiment ever being operated. Being at such shot dis¬ 
tance of the core lead to great difficulty in mitigating and 
controlling reactor induced accidental background. The 
high-accuracy measurement and subtraction of the huge 
gamma-ray reactor induced background was validated, 
as well as a novel method for assessing cosmic ray in¬ 
duced backgrounds at very shallow depth, correcting for 
the evolution of the atmospheric pressure during reactor 
OFF/ON periods. 

Eventually reactor antineutrinos at very short base¬ 
lines can be used to probe the Reactor Antineutrino 
Anomaly and to search for possible oscillations into ster¬ 
ile neutrino species. However, because of the current lack 
of precision in the emitted neutrino flux as well as the 
large accidental background level, no definitive conclu¬ 
sion could be obtained yet. A more refined evaluation of 
the expected neutrino rate would be necessary through 
further improvements of the detector response modelling, 
which are beyond the scope of this paper. Nevertheless 
Nucifer data can be now included into global analyses 
and it could perhaps constrain some of the still-allowed 
oscillation scenarios at short baselines. Dedicated experi¬ 
ments are clearly necessary to confirm or reject the short 
baseline oscillation hypothesis, such as the Stereo [36] . 
Neutrino-4 [37], PROSPECT [35], or SOX [33] 00] ex¬ 
periments expected in the next months or years. 
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The detection efficiency was improved by a factor of 3 
with respect to past reactor neutrino experiments dedi¬ 
cated to nuclear safeguards, reaching a mean daily rates 
of about 300l^e/day with a cubic-meter scale target vol¬ 
ume. Therefore, the achieved statistics is large enough 
to be of interest for cooperative monitoring regimes. Nu- 
cifer has also shown a great stability, suitable for relative 
monitoring of the neutrino rate within a few percent, and 
safe automatic operation for a few years with only little 
occasional maintenance. Hence, this operation regime 
approaches the requirements of any IAEA monitoring 
apparatus, making a Nucifer-like antineutrino detector 
suitable to monitor and safeguard nuclear reactors. 

As a first societal application, and within the frame¬ 
work of the Plutonium Management and Disposition 
Agreement for the disposal of Weapon-grade plutonium, 
we illustrated the possibility of monitoring the plutonium 
content in the Osiris core by showing that Nucifer could 


detect the presence of about 1.5 kg of Pu at 95% C.L, 
based on the actual data. This result would be improved 
with Nucifer deployed further away from a more powerful 
core, by intrinsically lowering the dominant gamma-ray 
accidental background component. 
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